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ABSTRACT 

As an alternative to the covariant Ostrogradski method, we show that higher- 
derivative relativistic Lagrangian field theories can be reduced to second differential- 
order by writing them directly as covariant two-derivative theories involving Lagrange 
multipliers and new fields. Notwithstanding the intrinsic non-covariance of the Dirac's 
procedure used to deal with the constraints, the Lorentz invariance is recovered at 
the end. We develope this new setting for a simple scalar model and then its ap- 
plications to generalized electrodynamics and higher-derivative gravity are outlined. 
This method is better suited than Ostrogradski 's for a generalization to 2n-derivative 
theories. 

PACS numbers: ll.10.Ef, ll.10.Lm, 04.60 
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1. Introduction 

Historically, field theories with higher order Lagrangians range from Higgs model 
regularizations [1] to generalized electrodynamics [2] [3] and higher-derivative (HD) 
gravity [4]. A procedure was later devised to reduce them, by a Legendre transfor- 
mation, to equivalent second-order theories [5] where a subsequent diagonalization 
explicitly displays the particle degrees of freedom [3] [6] [7] . 

The validity of the formal Lorentz covariant order-reducing method adopted there 
has been checked in an example of scalar HD theories by a rigorous study of the 
phase-space [7]. In this procedure, a generalization of the Ostrogradski formalism 
to continuous relativistic systems (2n-derivative because of Lorentz invariance) is 
carried out. In it, some of the field derivatives and the generalized conjugate momenta 
become, after a suitable diagonalization, new field coordinates describing the degrees 
of freedom (DOF) which were already identified in the particle propagators arising in 
the algebraic decomposition of the HD propagator. 

By using Lagrange multipliers, a variant of the Ostrogradski method for me- 
chanical discrete systems has been proposed which allows to show the quantum (Path 
Integral) equivalence between the modified action principle (first order Helmholtz La- 
grangian) and the starting HD theory [8]. For relativistic field systems, a similarly 
inspired procedure can be followed in which the multipliers let to write the HD theory 
from the outset as a second order (constrained) covariant one which lends itself to a 
particle interpretation after diagonalization. 

In this paper we implement this new setting by means of the use of Lagrange 
multipliers in a Lorentz invariant formulation of a relativistic scalar field theory. The 
Dirac method [9] prescribes the identification of the primary constraints arising in the 
definition of the momenta. These constraints are added to the starting Hamiltonian 
by means of new multipliers, and then they are required to be conserved by the time 
evolution driven by this enlarged "total Hamiltonian" through the Poisson Bracket 
(PB). This may give rise to secondary constraints, the conservation of which can in 
turn generate more secondary constraints. The process stops when constraints are 
obtained that can be solved by fixing the values of the new multipliers. We then 
use the remaining constraints to eliminate the starting multipliers and the momenta, 
ending up with a two-derivative theory depending on just its true DOF. Since the lat- 
ter appear mixed, a diagonalization works finally out the independently propagating 
DOF. 

As long as time evolution is analized, the true mechanical Hamiltonian (i.e. the 
energy) of the system must be used. Then one cannot benefit of the compactness 
of the Lorentz invariant procedures introduced in [7], so we are forced to deal with 
non covariant objects and face the diagonalization of larger matrices. The relativistic 
invariance of the system is recovered at the end of the process. 

From the methodological poin of view, the new treatment of HD scalar theories 
that we present here provides a sharp departure from the more traditional Ostro- 
gradski approach . Moreover, it is implementable and may prove advantageous in 
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generalized electrodynamics, HD Yang-Mills and linearized HD gravity as well. On 
the other hand, even if the explicit calculation already in the sacalar case becomes 
practically intractable beyond six-derivative order, we have partial results that lend 
themselves to generalization to arbitrary n better than the Ostrogradski method does. 

In Section 2 we treat n = 3, the case n = 2 being too much trivial for the 
illustrative purposes we pursue here. Some results regarding the extension to arbitrary 
n are presented in Section 3 which, whith the help of a plausible ansatz, allows the 
explicit calculation for n = 4. Section 4 discusses some possible applications of the 
approach to more relevant vector and tensor field theories. The Conclusions are in 
Section 5. An Appendix is devoted to a general inductive proof of the pure algebraic 
character (i.e. absence of space derivatives) of the secondary constraints. 
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2. n = 3 theory. 

We consider the six-derivative Lagrangian, 

C 6 = -^#1J[2][3J0- J0 , (2.1) 

where /i is an arbitrary mass parameter, [z] = (□ + mf) are KG operators, 
M = (12) (13) (23) , (ij) = mf — rrij > for i < j , and mass dimensions [/x] = 1 , 
[M] = 6 , [0] = 1 , [j] = 3 . 

As discussed in [7], (2.1) displays the general form of the free part of a higher- 
derivative scalar theory with nondegenerate masses m\ , , ms , the source 
term embodying the remaining self interactions and the couplings to other fields. 

C 6 can be reexpresed directly as a second-order theory with constraints, namely 

£6 = \ M Hal* + - + A 2 (^2 - [3]$,)] -jfo , (2.2) 

where ^3 = and Ai , A2 are Lagrange multipliers, so that (2.2) depends 
on five fields. Dropping total derivatives, in compact matrix notation, (2.2) reads 

£ 6 = -q T JC^ + -^ T M^> - J T $> , (2.3) 

where the vectors ^ and J , with components ipi , Ji , and the matrices 
/C and M. are 



V/x" 4 A 2 / 



so that [tpi] = 1 i = 1, . . . , 5 ; Ji = jS 3 i 
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M. is an operator with space derivatives present in Mf = mf — A 
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The canonical conjugate momenta are denned as 

^ = fr ■ (2-5) 
They are the components of a 5-vector IT for which one has 

n = /c^ . (2.6) 

Since /C is not invertible, not all the velocities ipi can be expressed in terms 
of the momenta and a primary constraint occurs, namely 

n 1 = ir 5 -ir 1 = , (2.7) 



as consequence of = ^0-^3 = 71 1 • There is only one such a constraint since the 



submatrix KL ab = 2M^ ab (a, 6 = 1, . . . , 4) is regular. In the following, indices 

a,b,... go from 1 to 4, while go from 1 to 5. The velocity ip§ is not 

worked out, and from (2.6) we have 

na = 2M >C ' ab ^ b + 2M 6a3 ^ a+2 ' ( 2-8 ' ) 
(do not sum over a), and therefore 

2M , 

tpa = —lC ab ir b - 5 al ip 5 . (2.9) 

The Hamiltonian is 



1 • • /i 6 • • 1 

H = ^a^a + 7T 5 ^5 - -1p a lC ab 1p b - T^V^ ~ ^iMij^j + J1p3 , (2.10) 

where %p a must be substituted according to (2.9). Then the dependence on ip 5 
cancels out and we have 

1 2M 1 

H = -^-^-^aKab^b ~ ^iMijlfjj + Jii>i . (2.11) 

In (2.11) only four momenta appear, together with the five fields ipi ; not all 
of the five momenta Tii are independent because of the primary constraint (2.7). 
The "total Hamiltonian", with five independent momenta, accounting for this is 

H T = H + (n 1 , (2.12) 
where ( is a Lagrange multiplier. 
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The stability of fii requires 



n l = {n 1 ,n T } PB = n 2 = ^((i3)i; 3 -^ 4 + ^ 2 ) = o . (2.13) 

This secondary constraint yields 

^4 = ^rV>3+^2 • (2.14) 

Further secondary constraints stem from the ensuing stability conditions 

O 2 = {O 2 ,H T } pB =O3 = (13)7ri-^ 2 7r 2 + ^ 2 7r 4 = , (2.15) 
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SO 



VT 4 = VT 2 - ^TT! , (2.16) 



and again 



n 3 = {n 3 ,n T } PB = n 4 = 



PB 

,6 



^(-(13) (23)^3 - ^ 2 (13>V>2 - A*Vi + A*Vs) = 



(once (2.14) has been used), from which one gets 



(2.17) 



, = <5f>,4^, . (2 , 8 ) 



The next constraint, after using (2.16), gives 



O 4 = {fi4,HT} PB =O5 = (13)(12)7r 1 -^ 2 (13)7r 2 -^ 4 7r 3 + 2^/C = , (2.19) 

and £ can be obtained as a function of 7Ti , 7T2 and 7r3 , thus bringing the 
generation of secondary constraints to and end. 

1~Lt being cuadratic in t/>'s and 7r's , guarantees an alternance of linear 
constraints involving the fields and the momenta. In spite of the occurrence of space 
derivatives in M , they cancel out and the constraints are algebraic. From this 
set of constraints, the multipliers ^4 and ^5 , together with their conjugate 
momenta 714 and tt^ , can be worked out. 
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The Lagrangian (2.3) can be expressed in terms of the independent variables 
ip a (a = 1,2, 3) . Notice that implementing these constraints in C 6 is legitimate 
as long as this operation does not erase the dependence of C 6 on the other variables. 
One obtains 

£ 6 = ip a IC ce pijj f 3 + ip a M a ptljf3 - jips , (2.20) 

where 



2M 
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(13) 




(13) 


(13)<23) 


A* 2 





(2.21) 



M a p = - (M a p + M aB M B p + NccaMap) 
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with a, /?,... = 1, 2, 3 ; A, £?,... = 4,5 ; and 



(2.22) 







1 



(13) 



(13) (13) (23) 

/2 ,,4 



■Va/3 = 

that allows to embody (2.14) and (2.18) in the closed form 



(2.23) 



4>A = J^Ap^p 



(2.24) 



The symmetric matrices /C and M. can be simultaneously diagonalized 
by the regular transformation 



where 



1pa = Ti a p(t>p 



(2.25) 



TlaP = 



(12)(13) 

M 4 
(13) 



r 




(2.26) 



The 3rd. row of the non-orthogonal matrix 1Z in (2.26), has been chosen so 
as to yield the correct source term in (2.29) (see below). The remaining six elements 
are uniquely determined by requiring 1Z to diagonalize IC and M. 
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The diagonalized matrices are 

tt.^^M^ffl,-^) , (2.28) 
so that (2.20) finally writes 

This shows again the Ostrogradski-based result [7] of the equivalence between the 
second-order theory (2.1) and the LD version (2.29) that reproduces the propagator 
structure. 
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3. Theories with arbitrary n. 

The general Lagrangian 

C 2n = -I^[l][2]... , (3.1) 

where M = Y[ (ij) , and (3 = n(n — 3) + 2 for dimensional convenience, can be 

i<j 

dealt with along similar lines. The 2-derivative constrained recasting of (3.1) is 

^ 2n = ^h^I 1 ^l + Al(^l-[2]^2) + ... + A n _ 1 (^n-l-M^n)]-^n , (3.2) 

with ip n = (f) , and Ai, . . . , A n _i being Lagrange's multipliers. In order to have 
a more compact notation we define 

ipa = /U" 2(n_Q) ^a a = l,...,n 

(3.3) 

ipA = A 4 2a A a A = n + a; a = 1, . . . ,n— 1 
so that [V>i] = 1 (i = 1, ■ • • , 2n - 1) . Then 

£ 6 = 1^ T /C^ + i^ T .M^ - J T ^ (3.4) 

with Ji = jSin , and the (2n — 1) x (2n — 1) matrices /C and are given 

by 

^Qj = &(5i,j-n+l + <5j',i-n+l) 

(3 5) 

M i:j = a [-(Mf5 hJ - n+1 + Mj5 Jyl _ n+1 ) + ^ 2 (5 hJ . n + <^_ n )] , 

n(n— 1) 

with a = M 2M — . This "mass" matrix contains again space derivatives. Here, 
as before in (2.8) and in the following, an underlined index means that Einstein 
summation convention does not apply. The canonical conjugate momenta are now 

dC 2n 

7Ti = —y- , (3.6) 

n = KM . (3.7) 



i.e., in closed notation, 



Defining the matrix /C 



< 6 = -/C a6 (a,6=l,...,2n-2) , (3.8) 



a 



one sees that det/C ^ , while det/C = . This means that we only have one 
primary constraint, namely 

fil = 7V 2 n-l ~ TTl = . (3.9) 
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Then tp2n-i is n °t worked out, while ip a (a = 1, . . . , 2n — 2) can be expressed 
in terms of 7r a and ip 2n -i • The first 2n — 2 components of eq. (3. 7), namely 



7r a = aJC ab tp b + aS an ipa+n-i (3.10) 

give 

^ a = -JC ab n b - 5 al ip 2n -i ■ (3.11) 
<j 

After checking that the terms in ip 2n -i cancel out, the Hamiltonian has the 
simple expression 

W = ^a7r a /C a6 7r b - ^ipiMijifjj + jip n . (3.12) 

In TL only 2n — 2 momenta 7r a occur against 2n — 1 fields ipi , be- 
cause of the primary constraint (3.9). One may restore the dependence on 2n — 1 
momenta by introducing the "total Hamiltonian" 

H T = H + (n 1 , (3.13) 

where ( is a Lagrange multiplier. 

From the stability condition on fii , a cascade of secondary constraints follows, 
eventually ending with an equation that determines the value of ( . We outline 
here the steps closely following the lines of section 2. 



Q 1 = {Q 1 ,'Ht} pb = ^2 = => V'n+l = V'n-l + ■ (3.14) 

Then 



n 2 = {tt2,H T } PB = n 3 = ij J \ 2n _ 2 + (ln}Tr 1 -Li 2 TT 2 (l-5 n2 )-2a(;5 n2 = . (3.15) 

If n = 2 , eq.(3.15) gives ( in terms of tti and ir 2 , and the cascade stops 
here. If n > 2 , eq.(3.15) yields 

VT 2 n-2 = 5-7T1 + 7T 2 • (3.16) 

LI 2 

The next step is fi 3 = {fi3,7ir} ps = fi 4 = , which together with (3.14) 
gives 



tP n+2 =^ n _ 2 + ^((ln) + (2,n-l))Vn-i + ^(ln)(2n)^ n , (3.17) 
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and, proceeding further, we obtain for the momenta 



4 = {^4,'Ht} pb =^5 = ^ 4 ?T2n-3 ~ (In) (1,71 - 1)tTi + 

+ fi 2 ((ln} + (2,n-l})n 2 -fi 4 7T 3 (l-5 n3 )-2afi 4 C5 n3 = , 

where (3.16) has been taken into account. Again, if n = 3 , the process stops here 
and we have reproduced the results of section 2. If n > 3 , eq.(3.18) yields 

7i2n-3 = -^(ln){l, n - l>7ri - \{{ln){2, n - l))vr 2 + tt 3 , (3.19) 

and the process goes on. 

For ilustrative purposes, we complete here the steps that cover the case n = 4. 
fi 5 = {Q 5 ,n T } PB = ^6 = , yields 



t/v+s = ^ n _3 + «3, n - 2) + (2, n - 1) + (ln))^„_ 2 + 
A* 

+ -1 ((2, n - 1)(3, n - 1) + (ln)((2, n - 1) + (3n)))V„-i+ (3.20) 
A* 

+ -l(ln)(2n)(3n)V„ , 
A* 

and 06 = {^6)^t} pb = O7 = gives £ in terms of 7Ti , 7r2 , ^3 and H4 . 

In general, for a fixed n , the quadratic dependence of H on 7r$ and 
, together with the primary constraint Oi , leads to a set of secondary con- 
straints fifc that splits in two classes according to k being even or odd. A con- 
straint Q,2j (j = l,...,n — 1) is linear in ipi and gives ip n +j i n terms of 
ip n , . . . , ift n -j. A constraint &2j-i (j ' = 2, . . . , n — 1) is a linear combination of 
7Ti and gives 7T2 n -j in terms of 7Ti,...,7Tj . Finally, f^n-i fixes the value 
of C and stops the process. 

One can prove that the constraints ^2n-i do not contain space derivatives, 
even though the elements of M. involved in their calculation contain the Laplacian 
operator. This will be shown in the Appendix. 

Like in (2.24), we take 

V>a = Ma^h (3.21) 
with indices (a, /?,... = 1, n) and (A, B, . . . = n + 1, . . . , 2n — 1) 



11 



Now Af is a (n-l)xn numerical matrix whose three first rows can be read 
from (3.14), (3.17) and (3.20). One sees that the elements of the j — th row became 
more and more complicated for bigger j ; we have always in the n — j — 1 
first places of this row and 1 in the n — j position (see Appendix), although we 
lack a general expression in closed form except for some elements. 

Then, the Lagrangian again is 



C 2n = + 1p a M af 3ll)f3 ~ j^n ■ (3.22) 

The n x n matrices K, and M. have the same structure in terms of /C , 
M. and Af given in (2.21) and (2.22). However, the difficulty of explicitely finding 
the elements of Af , makes IC and M hard to calculate. 

The diagonalization of (3.22) will be accomplished, as in (2.25), by a n x n 
real matrix 1Z . We again impose TZ n p = 1 {(5 = 1, . . . , n) to ensure the cor- 
rect form of the final source term. The requirement of simultaneously diagonalizing 
K, and A4 , yield n(n — 1) quadratic equations that determine the n(n — 1) 
remaining elements of 1Z . The existence of such a regular 1Z with real elements 
is by no means guaranteed a priori, but the results of [7] , showing the equivalence of 
the HD and the LD Lagrangians, spurs us to look for it. Although obtaining 1Z is 
almost impossible already for n = 4 , we guess his general form, namely: 



lZ a/3 = 1 ; (a = n) 

^ = (-l)->- 2 ^- a )(/3,a + l)(/5,a + 2)...(/?,n) ; {(3 < a < n) (3.23) 
1Z aP = ; (a</3) . 

Of course, this TZ is just (2.26) for n = 3 . For n = 4 in fact, from 
(3.14), (3.17) and (3.20) the matrix Af is known, and assuming (3.23) we obtain 
the LD Lagrangian 



2 (1) 2 (2) 2 (3) 2 (4) ( 3 _ 2 4) 

-j(<f>l +02 +<f>3 +04) , 

where (i) = j^M Y[ j^yy , and M is given in (3.1). This is the result expected 

from the covariant Ostrogradski method shown in [7] . This success strongly backs the 
ansatz (3.23). 

Finally, we want to remark that the case n = 2 is trivially contained in the 
general n case considered in this section. 
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4. Applications to other theories. 



The constraint method we have developed for scalar theories can be implemented 
for HD vector and tensor theories as well. 

In the case of HD vector theories, a most general example is the gauge-fixed 
generalized QED, given by 



C — —-F F^ v - F nP" 
L ~ 4m 2 ^° 



(4-1) 



The structure of its constraints has been studied in [3] by a canonical forcefully non- 
covariant analysis carried out on the Ostrogradski-based order-reduction procedure. 

A recasting of a higher-derivative gauge-invariant Yang-Mills theory as a two- 
derivative one by means of constraints has been done in a non-covariant 3+1 way 
[10], while we are interested in keeping the explicit Lorentz covariance at this stage. 

Dropping total derivatives, (4.1) may be written as 

L = \ A "K + C 2 <)n(n(^ + ^) + V P »)A»-j»A» (4.2) 

d d d d 

where = r]^ — - J ^ L and uj^ v = - J ^ L , with Minkowski metric 

Tj^ v = diag(l, — 1, — 1, —1). Then, omitting indices, the covariant two-derivative con- 
strained version may be readily written as 

C = \{A □ (9 + eu)B + A(B-(n(^ + ^) + v)A)}-jA (4.3) 

where the four-vector A is the multiplier and B is a new vector field. The Lagrangian 
(4.3) is local and is regular in the time derivatives of the fields involved. Therefore 
it is adequate for defining conjugate momenta 7r^, 7r^ and tt 1 ^ upon which the Dirac 
method and subsequent diagonalization can be implemented. 



The covariant Ostrogradski order-reduction of the four-derivative gravity leads 
to a two-derivative equivalent in which the particle DOF can be fully diagonalized [6]. 

The constraint technique for the order-reduction of a pure four-derivative confor- 
mally invariant gravitational Lagrangian has been already used in a 3+1 non-covariant 
form [11], where further first class constraints from Diff-invariance occur. In a covari- 
ant treatment and for the general case including also two-derivative terms [12], a 
seemingly similar method is adopted where in place of the Lagrange multiplier a less 
trivial auxiliary field featuring a squared (mass)term is used. A little work shows 
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however that this method is identical to the covariant Ostrogradski's [7], as shall be 
discussed elsewhere. We illustrate this here on the grounds of the scalar model. 

Consider the four-derivative Lagrangian, 

>C 4 = -^0[1][2]0-J0 , (4.4) 

which also reads 

>C 4 = -^[p0 2 + S 0(n0) + (D0) 2 ]-i0 , (4-5) 

where p = m 2 m 2 and s = m\ + m 2 . The covariant Ostrogradski method, in a 
sligthly less refined version that the one presented in [7], would define a conjugate 
generalized momentum tv = • The Legendre transformation performed on it 

leads to a Hamiltonian-like density from which the following two-derivative Helmholtz 
Lagrangian is derived 

C H = txU4> + 1(12) 7T 2 + 1(12) (P 2 + 7T0 . (4.6) 

On the other hand, by using the auxiliary field technique of [12] the higher-derivative 
term is brought to second order by writing (4.5) as 

/: 4 = -l^y[p0 2 + ^(n^) + A(n^)-lA 2 ]-^ , (4.7) 

where the equation of motion for A recovers (4.5) when substituted back in (4.7). 
Now, in spite of their quite different look, (4.6) and (4.7) are related by the simple 
field redefinition tv = — ^ -^j + s 4>) ■ 

The covariant constraint method introduced in this paper provides a new ap- 
proach. The most immediate application in higher-derivative gravity regards the 
linearized theory, usually considered when analizing the DOF. Take for example the 
four-derivative Lagrangian 

C = V=g[aR + bR 2 + cR liV R fU/ ] . (4.8) 

The linearization around the flat Minkowski metric, namely = r/^ + h^ v , 

simplifies it to 

C-=\h ilv [{aG + 2bfiUP + 2cf 1 UP)UP}^ pa h p(J , (4.9) 
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where f}^' pa = \{ff f, t) va + T}^T] vp ) , ff v ' pa = v ^ri pa , G = \f\ - f) and 
Pfiu, P a = liw^pVva+Uvprj^-CLi^rjpcr-fj^pv) . Omitting indices, the order-reduction 
of the theory by means of a Lagrange multiplier yields the two-derivative local La- 
grangian 

C=^[h(aG + 2bf l nP + 2cf ] nP)f + A(f-nPh)} , (4.10) 

where is a new field and A Mi , is the multiplier. Of course, because of the 

Diff-invariance, first class constraints will remain when the Dirac procedure is carried 
out. 

5. Conclusions 

We have shown how to deal with 2n-derivative relativistic scalar theories by 
writing them directly as second-order constrained Lagrangians with more fields and 
suitable Lagrange multipliers. The corresponding canonical conjugate momenta are 
subject to primary constraints, whose conservation in time gives rise to a finite chain 
of secondary constraints according to the Dirac's procedure. Though expected, a non 
trivial result is that these constraints, later used to extract the final DOF, are purely 
algebraic relations that do not involve the space derivatives. 

Once the constraints have been implemented, we are left with a second-order 
Lagrangian for the DOF of the system. We have performed explicitely the diagonal- 
ization for n=3 and used an ansatz to work out the case n=4, confirming the result 
obtained in [7]. Though not proven , this ansatz is given a plausible expression for 
arbitrary n, namely (3.23). This step towards the explicit generalization to higher n, 
gives this method an advantage over Ostrogradski's. 

The applications to more interesting theories like HD generalized electrodynamics 
and HD Diff-invariant gravity illustrate also the fact that the order-reducing methods 
used in the literature fall in two categories: the one based in the covariant Ostrogradski 
and the one based in the contraints by Lagrange multipliers. The methods based on 
auxiliary fields with a quadratic term, which may look like a variant of the multipliers, 
actually belong to the first category and have no obvious extension beyond the four- 
derivative order. 

In vector and tensor field theories where gauge symmetries occur, the correspond- 
ing first class constraints live together with the second class ones worked out in this 
paper and survive the order-reducing procedure as long as gauge fixings are not con- 
sidered. The method may then prove useful for a detailed analysis of the constraints 
from gauge (or Diff-)invariance in these HD theories, chiefly of the fate of the scalar 
and vector constraints of Hamiltonian gravity. 
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Appendix 



We prove, by induction, that the constraints VI23 (j = 1, . . . ,n — l) involving 
the fields, do not contain space derivatives because the Laplacian operators cancel out. 

One first sees, by inspection, that this statement is true for O2 : the calcula- 
tion leading to (3.14) is 

fi 2 = a [n 2 ^n-i + {Ml - Ml)^ n - /i 2 ^ n+ i] = , 

(A.l) 

where the cancelation of the Laplacian operator is apparent, i.e. 

M 2 -M 2 = m 2 -m 2 = (ln> , 

(A.2) 

and obviously no summation is understood for repeated indices. Then, let us suppose 
that, after taking into account the preceding constraints, one has that in the constraint 

^2a = a[^ 2a ^ n _ Q +ai^ n _ Q+ i+a 2 ^n-a+2 + - • -+a a -i ip n -i +a a ip n - fj? a ip n+a ] = , 

(A.3) 

for a = 1, . . . , n — 2 , the coefficients a±, . . . , a a are real numbers, as are those 
found in (3.14), (3.17) and (3.20). We now prove that this is also true for 2q +i . 
In fact 

^2a+l = /i a 7T2 n -a-l + 0>1^2n-a + CL2^2n-a+l + ■ ■ ■ + da-1^2n-2 + a a 7Vi — 

- [i 2 7l a+ l(l - <S n ,a + l) - 2(7^ 2 C<5n,a+l , 

(A.4) 

from which 

^2(a+l) = ^2a+l = {^2a+l > } pB = 

= <7[- / U 2Q M 2 _ a V'n-a - aiM 2 _ Q+1 -i/>n-a+l - ... - a a _i M 2 _ t 1p n _ l - 
- (l a M 2 J) n + l? Q : M 2 +a lP n+a + ^ 2 ( Q + 1 Vn-a-l + ^ai^n-a + • • • + 

+ ^ 2 a a _i^ n _2 + n 2 a a ip n+1 - ^ 2(a+1 Vn+a+i] = . 

(A.5) 
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The crucial point now is that, when working V'n+a out of (A. 3) and substi- 
tuting it in (A. 5), only differences of squared masses Mf occur as in (A. 2), thus 
cancelling out the operators A . Then, by substituting also ip n +i from (A.l), 
one gets i/j n+a+ i as a sum of linear terms in ip n , ip n -i, . . . , VVi-a-i > the coeffi- 
cient for the last one being the unity. This ends the inductive proof. 
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